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Qh' Abstract 

I The influence of disorder on the transmission through periodic waveguides is 

^2 ' studied. Using a canonical form of the transfer matrix we investigate dependence 

of the Lyapunov exponent 7 on the frequency v and magnitude of the disorder 
a. It is shown that in the bulk of the bands 7 ~ o"^, while near the band edges it 
has the order 7 ~ a^l"^ . This dependence is illustrated by numerical simulations. 

> 
(N 

1 Introduction 

O I Propagation of time-harmonic electromagnetic waves in one- dimensional periodic op- 

tical waveguides can often be reduced to the eigenvalue problem 



■ n^[x) dx 

■ where iplx) is a Cartesian component of electric or magnetic field, n{x) = n{x + 
i) is the index of refraction which we suppose to be periodic. A way to find the 
frequencies u for which waves can propagate is by computing the so-called characteristic 
or transfer matrix M that maps two-dimensional vector with components ip and ip' at 
the beginning of the period to the value of that vector at the end of the period 



M 



(2) 



Matrix M is an analytic function of the frequency u and detM = 1. Propagation 
frequencies u of the problem (1) satisfy the condition 

|trM|<2. (3) 

If (3) does not hold for some frequency u then the corresponding wave cannot propagate. 
The spectrum of operator (1) (which is the closure of the set {i^^} for all propagating 
frequencies) has a bandgap structure with the bands defined by |triVf| ^ 2. 
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In reality, however, periodic waveguides contain finite number of periods A^. Equa- 
tion fl) in tliis case becomes 



(Pi) 



(4) 



wfiere 



Unix 



(5) 



1, x^[0,Ni], 
n(x), X e [0,Ni]. 

Sucli waveguide supports any frequency z/ ^ 0, but transmission for different frequencies 
varies. The transmission and reflection coefficients tN{v) and r7v(z/) are determined by 
solving the scattering problem 



dx"^ 



+ v n^{x)'ilj{x) = 



with 



ifji^x) 



e*'^^- + e-*'^^ x < 0, 



X > Ni. 



(6) 



(7) 



The transfer matrix over the interval (0, Ni) is M{Nt) = M'^ . The transfer matrix 
Mp in the Priiffer form maps two-dimensional vector with components ip and v at 
the beginning of the interval to the value of that vector at the end of the interval, i.e. 



Mp{Ni) = T~^M^T, T 



\i^\ 



-1/2 









.|l/2 



This matrix defines the transmission coefficient [1,2] 

4 4 



t 



N\ 



\M^{Ni)f^^ + 2 \\T-^M^T\ 



where ||iVf ||hs is the Hilbert-Schmidt norm of the matrix M = [Mij] 



\M\ 



tr MM* = ^M, 
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Random imperfections in periodic waveguides can significantly affect their proper- 
ties, especially near degenerate band edge, see [3, 4]. This paper deals with the influence 
of disorder on the propagation properties of the waveguide. Let M^, k = 1,2, . . . , N 
be transfer matrices of the disordered waveguide on the periods [{k — l)i,ki]. Then 
formula (8) becomes 



M = ^ 

T-^MnMn^i ■ ■ ■ MiT 

The norm of the product of random matrices is related to the Lyapunov exponent 

ln||M^Mjv-i...Mi||Hs 



(9) 



7 = lim 

Af-s>oo 



N 



(10) 
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If this limit exists then 



\tN\ 



-■jN 



N ^ oo. 



fill 



Thus, the Lyapunov exponent allows one to estimate the energy transmitted through 
N periods of the waveguide. In the paper we study 7 as a function of disorder strength 
a and show that it has different behavior depending on the location of frequency u in 
the band. 

We will consider two types of disorder. Disorder of the first type concerns pertur- 
bation of physical properties of a periodic waveguide and is described by the equation 



dx"^ 



0, A{x + i)= A{x), B{x + l) = B{x), (12) 



where is a stationary random noise with correlation length A ^ 1, E^a 
|.^a| < 1, and a <^ VA, i.e. converges in law to the white noise as A — )■ 



Disorder of the second type is related to perturbation of the length of the period. 

We shall show that in both cases matrices Mk can be represented (c.f. [5]) as a 
multiplicative perturbation of M 



Mk = M{l + aVk + 0{a^)) 



where tr = 0: 



Ck —C,k 



(13) 



(14) 



and entries C,k,VkXk of 14 in the limit as A — )■ form Gaussian random vectors 
independent for different k. The covariance matrix 



^k ^kVk ^kCk 



^kVk 



VkCk 



^kCk VkCk Ck 



(15) 



of elements Ck,Vk,Ck is non-degenerate for small A > 0; moreover, the limit Bo = 
limA-j.o-BA exists and det-Bo > 0. Later we will refer mostly to (13) rather than a 
specific form of the equation. Thus our results can be applied to more general periodic 
media with small and short correlated disorder. 

The paper contains the following results. Section 2 concerns the deterministic 
problem. We discuss bandgap structure of the spectrum of the periodic problem (1) 
and properties of the transfer matrix M near the fixed band edge. We provide an 
example showing that unlike the Schrodinger equation, the gaps of the optical equation 
(1) are located irregularly and their length does not go to zero. After that throughout 
the paper we will consider u only in a fixed band. Then we show that the transfer 
matrix in a neighborhood of the band edge can be reduced to a simple canonical 
form defined by equation (1) with n = 1. In Section 3 we show that the canonical 
representation of the transfer matrix holds for two major cases: the disorder of the 
refraction index and length of the period. Then in Section 4 we estimate the Lyapunov 
exponent of randomly perturbed periodic waveguide by the Lyapunov exponent of 
the simplest model - the Schrodinger equation with white noise potential. Section 5 
contains asymptotic analysis of the white noise model for the Schrodinger equation 
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and derivation of expression for the Lyapunov exponent. We show that 7 ~ cr^ in the 
bulk of the band and 7 ~ o^l"^ near non-degenerate band edges, and present resuhs 
of numerical simulations. This asymptotics for the Schrodinger equation leads to a 
similar estimate for the optical problem due to results of Section 4. 



2 Spectral properties of periodic optical waveguides 

We start with an example showing the main difference between the optical and Schrodinger 
equations. Let the index of refraction n{x) be ^-periodic, piecewise constant with two 
different values ni,n2 on two successive subintervals of length li,l2 on the interval of 
periodicity i = li + 12- Then 

M = M2M1, (16) 



where 



Mi 



Hence 



, , , sin(?7,j/jZ/) 
cos(nj/ji/) 

-riiU sin{niliu) cos{niliu) 



1,2. 



(17) 



tr M = 2 cos(?T,iZiz/) cos{n2l2i^) 

1 smlniliu) sm(n2<2Z^) 

n2 nij 

A cos [av] — B cos [hv) , 



18) 



^ A ^ f ni n2\ , „ ^ f ni n2\ . , , , , , tp 

with A = -\ 1 +1, B = -\ \ -1, a = ni/i + 71.2/2, = Uili - n2L2- it 

2 \n2 nij 2 \n2 rii J 

numbers rzi/i and 712/2 are rationally independent then almost periodic function /(z/) = 
tr M oscillates along the u axis between A + B = — -\ — - > 2 and —(A + B) < —2. 

7T,2 111 

In this case there are infinitely many gaps in the spectrum and their length does not 
vanish as z/ — 00. In fact, there are values z/ = z//^. — > 00 as i — t- 00, such that 
cosaz/fc = 1, cosbuk ~ —1, f{i^k) > \{A + B + 2) > 2 and fu is bounded uniformly in 
z/. At the same time the gaps can be arbitrary small, see Fig. 1. 
We remind that in the case of periodic Schrodinger equation 



+ V{x)i){x) = v'^i){x), V{x) = V{x + i), 



(19) 



the gaps are located periodically (see Fig. 2) and asymptotically vanishing [6]. For the 
7T,-th gap Qn one has 



e„, 5„, — > as 71 00. 



(20) 



Asymptotic behavior of en,Sn, depends on the smoothness of potential V{x) (see [6]). 

Our next result concerns the structure of the transfer matrix near the band edges. 
Let us consider one of the spectral bands /3 = [z/q, ui] and analyze the transfer matrix 
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Figure 1: Graph of the discriminant curve of the periodic optical waveguide with 
li = 2, I2 = 0.2, ni = 77-2 = 6 oscillates between A + B = — -\ — - and —{A + B) 

77-2 rii 

exhibiting infinite number of bands and gaps as z/ — )■ 00. 



M as a function of u when u ^ uq. The case z/ — ?■ z/i can be considered similarly. 
Inside the band /3 we have 

^trM{u) < 1. (211 



The eigenvalues of Af have the form e 



< 1. 

= ulu), and hence 



cos a; = - tr M{iy). 



(22) 



This matrix can be reduced to the diagonal form for any 1/ strictly inside the band. At 
the band edges, however, both eigenvalues of M equal either +1 or —1. Throughout the 
paper we consider only non-degenerated band edges, i.e., we assume that the derivative 
of the function tr iVi"(z/) is not zero at the band edge (where |trAi"(z/)| = 2). It means 
that we do not consider points like u = 7.755 in Fig. 1. Equation (22) immediately 
implies that u is an analytic function of ^/u — uq in a neighborhood of a non-degenerate 
band edge i> = z/q. It is also important that the transfer matrix cannot be diagonalized 
at a non-degenerate band edge, but rather reduced to the Jordan form. 

The simplest example of the discussed situation is given by the equation ip" + uj'^ip = 
0, ^ X ^ 1 . The transfer matrix for that equation is 



M(u) 



cos a; 

-oj sin u cos oj 



(23) 



and it turns into the Jordan block at the non-degenerate band edge u = 0. The next 
Proposition shows that arbitrary transfer matrix in a neighborhood of a non-degenerate 
band edge can be reduced to that simple form uniformly in u. 
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Figure 2: Graph of the discriminant curve of the periodic Schrodinger equation (sohd 
hne) with li = 2, I2 = 0.2, rii = 4, n2 = 9. Dashed hne corresponds to 2cosi^£. 



Proposition 1. Let M = M{u) be a transfer matrix of equation (1) in a neighborhood 
of a non- degenerate band edge v = uq, and let e^^'^ be the eigenvalues of Mlu), u = 
u}{i'), a;(i/o) = 0. Then there exists an analytic matrix-valued function D{uj), detD ^ 
0; such that 

' cosw ^ 



M = D-^M(iy)D 



-u sm u cos u 



(24) 



Remark. This proposition is vahd for general matrices with the same properties as 
those mentioned above for Miu). 



Proof. Let 



M(u) 



a b 
c d 



2 ^ ' 



a -\- d 



cosw, 



(25) 



where u = a;(z/) is real for u > uo (in the band) and u = w(z^) is complex for u < uq 
(in the gap). We have that \b\ + |c| > for u = uq, since M^u) is not diagonalizable. 
We may assume that b{uo) 7^ 0, and therefore \b{h')\ ^ 6 > for |z/ — z/o| ^ e. 

The eigenvectors of matrix M corresponding to the eigenvalues e^*'^ have the form 



'ViM 



a — e 



Vectors Vi 2{iy) form a basis for any u > 0. However, Vi = V2 when u 
basis degenerates. 

In order to avoid the degeneracy, we introduce new basis 



(26) 



0, and the 



Ul = -[Vi+V2), U2 



2iuj 



.Vi - V2 . 



(27) 
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The new basis is uniformly non-degenerate, and matrix M in this basis has the form 

(28) 



i.e. (24) holds with 



smcj 



COS 00 

UJ 

-u sin u cos u 



D 



a — cos u 





sincj 



(29) 
□ 



3 Transfer matrix for disordered wave guide 

Here we will show that our main assumption (13) on the form of the transfer matrix 
holds for typical random perturbations of deterministic equations. Two classes of 
problems will be considered. The first one concerns the equation 



+ [z/M(x) + B{x)] ip = 0, A{x + £)= A{x), B{x + = B{x), 



and its random perturbation 



dx'^ 



+ [u^A{x) + B{x)] ij + a^4^^lj = 0. 



(30) 



(31) 



Here is a stationary random process with short correlation length A <^ E^a = 0, 
which converges to the white noise as A — 0, and a <^ -\/A is a small parameter 
(strength of disorder). 

Let ipi{x), ip2{.x) be the fundamental solutions of equation (30) with ipi{Q) = 
1, 'ip'iiO) = 0, 4'2{0) = 0, 4'2{0) = 1, and let ipii^x), 4'2{x) be similar fundamental 
solutions of the perturbed equation (31). The transfer matrices on the interval [0,i] 
for (30), (31) are given by 



M(z/) 



Miu) 



(32) 



Transfer matrix Nl^iv) on the interval [/c£, ik^ 1)£] is defined similarly (matrix Miv) 
does not depend on /c). 

Lemma 1. The 'perturbed transfer matrix Mk{i') of equation (31) on the interval 
[k£, {k + 1)£] has the form 



M, = M{l + aVk + Oia^)). 



(33) 



Here matrix Vk has the form (14), trV^' = 0, and the covariance matrix (15) has 
a limit as A — > 0, = B + 0(A), where B (and therefore B^ for small A) is 
non- degenerate. 
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Proof. Without loss of generality we assume that k = and omit subindex k in the 
proof. We compute iVf(i/) using first order perturbation theory. Solutions ipi^2 for the 
disordered waveguide can be found from the integral equation 



i/ji{x) = ipi{x) + a I G{x, s) ^^jlJ ilJi{s) ds 



where 

1pl{s)ij2{x) - iJi{x)ij2{s), S <X] 

0, s > X. 

Solution of (34) in the linear approximation has the form 



G(x, s) 



tpi{x) =ipi{x) + a / G{x, s) 



iJi{s) ds. 



If we put 



ipi{s)ipj{s) 



Us) 



ds. 



1,2, 



(34) 



(35) 



(36) 



(37) 



then the transfer matrix of the disordered waveguide becomes 



M 



m) ^2 



vum) - vi2m)) vum) - V22m) 



+ 0(a' 



M + CrM 



-Vl2 -V22 
VU Vl2 



O(a^) = M{I + <jV + 0{a'^)) 



-Vi2 -V22 



(38) 



Random 



where M is the unperturbed transfer matrix and V 

rill Tli2 

variables rjij are asymptotically Gaussian when A —t- as it follow from the central 
limit theorem. Their covariance matrix B/\ (15) has a limit as A — 0, B/^ = B+0{A), 
where B is expressed in terms of the fundamental solution as follows 



B 



lo i'I{s)i!l{s) ds iJi{s)ijl{s) ds - ij'f{s)ij2{s) ds 

ipi{s)ipl{s) ds Io'^t{s)ds - ilJi{s)ipl{s) ds 

■ ipf{s)ip2{s) ds - ilJi{s)ipl{s) ds Io'^f{s)ds _ 



(39) 



Observe that B is the Gram matrix for the system of functions —tpKs), ipi{s)ip2{s) , 
■j/^Ks) on [0,i]. The degeneracy of B would imply that this system of functions is 
linearly dependent, i.e., 



cniplis) + Ci2Ms)Ms) + C22ip2is) = 0, se [0,i], 



(40) 



for appropriate constants Cii,Ci2,C22- The latter is possible only if ipi and ip2 are pro- 
portional, which contradicts to linear independence of ipi and ?/'2- Thus, the correlation 
matrix B is non-degenerate and so does Bj\ for sufficiently small A. □ 
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Let us show that (13)-(15) holds for another important class of problems related to 
the wave propagation in a layered medium with layer thickness disorder. Consider a 
periodic waveguide described by (1) with ^-periodic piecewise-constant function n{x). 
Each period [ki, {k + l)i] consists of m subintervals of the length /j, £ = /1 + /2 + . . ■ + lm, 
where n{x) is equal to rii on i-th subinterval. Suppose that each subinterval changes 
randomly to + (T^j), where < 1 are independent random variables with densities 
(for example, C,i are uniformly distributed on the interval (—1,1)) and cr is a small 
parameter. We assume that distributions of C,i do not depend on k and omit index k 
below. Then the transfer matrix of z-th subinterval becomes 



M, = 

where Mi 
the form 



cos riiuli^l + a^i) 



sinniZ/Zi(l + aE^i 



riiU 



—•riiV sin niiyli(l + a^i) cosnji^/j(l + a^i) 



M,iI + aV, + Oia^)), (41) 



Ml 



(7=0 



is the unperturbed transfer matrix of i-th subinterval, and Vi has 



1 




trK = 0, 



1, 2 m. 



(42) 



(45) 



The transfer matrix for the period (consider for simplicity the case m = 3) is 
M = M3M2M1 = Ms{I + aVs)M2{I + aV2)Mi{I + aVi) + 0{a^) 

= MgMaMi [l + a {Vi + Mf V2M1 + Mf ^Ms" V3M2M1) + 0{a^)] . (43) 

The last two terms in the parentheses are similar to the matrices V2 and V3, respectively, 
and therefore by (42) trVi = trMfVsMi = tr M^^M^^V3M2Mi = 0. Hence, 

M = M {I + aV + 0{a^)) , (44) 

where trV = 0: 

c 

A similar result is valid for arbitrary number m of subintervals on a single period. The 
entries ^, 77, ( are a linear combination of random variables ^j, i = 1, 2..., m. 

Let us show that the covariance matrix B of the vector p = [C,, t], Q is non- 
degenerate. Each matrix Vi contains only one random variable C,i, and p has the 
form p = Yl^i ^i^ii where are three-dimensional vectors defined by nonrandom 
matrices. 

Without loss of generality we may assume that there are at least three linearly 
independent vectors among {aj}. If m ^ 3 one needs only to note that two adjacent 
subintervals have different refraction indices rij 7^ n^+i. Also, we double the period if 
m = 2. Now 

(m m \ m 

aiii ® aj^j = E (e- ) Y'^i® ^ 7^ 0- (46) 

i=l j=l J 1=1 

We cannot expect that the law of p is close to the Gaussian. However, if we choose 
N£ as a new period with ^ 1 and assume that Na^ ^ 1 then the distribution of 
^i^i ^^^^ tie close to the Gaussian. Thus, an analog of Lemma 1 holds for the 
case of layer thickness disorder. 
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Comparison theorems for the Lyapunov expo- 
nents 



The goal of this section is to compare the Lyapunov exponents 7 for the optical waveg- 
uide (31) and for the Schrodinger equation with the white noise potential. The latter 
equation can be obtained from (31) with A = —1, B = when A — ?■ 0, and it has the 
form 

-ij"{x) + aw{x)ij{x) = uj'^ij{x). (47) 

Here we use different notation for the spectral parameter in order to compare these 
two problems with 00 = 00(1'). 

Proposition 2. Lyapunov exponent 71 = 71(0", i^) for optical equation (31) with v 
in a neighborhood of a band edge uq and Lyapunov exponent 72 = 72(0", w) for the 
Schrodinger equation (47) with u = (X'(z/) defined in (22) and (24) estimate each other 
uniformly with respect to a and v, i.e. there exist constants Ci, C2 such that C172 ^ 71 ^ 
C2I2 for all < a < (Tq, |z/ — z/qI < /iq for some ctq, /io- 

The transfer matrices Mk for the optical problem have the form (13), where matrix 
M is similar to M in (24): 



cos a; 

-00 sin 00 cos 00 



{l + aD-'VkD + 0{a^)) , uj = uj{iy), (48) 



with Vk satisfying (14), (15). 

The transfer matrices iVffc(z/) for the Schrodinger equation (47) also satisfy (13)-(15) 
with M = M, since A = -1 and B = 0: 



cosw 
—00 sinu 



cos a; 



{l + aVk + 0{a^)). 



(49) 



Here Vk are different from those in (48), but both satisfy (14), (15). Now, to justify the 
proposition one needs only to note that the similarity transformation of all matrices 
Mk (with the same matrix D) does not change 7 (see section 2) and apply the following 

Lemma 2. Consider two sequences of independent identically distributed random Gaus- 
sian matrices 



Mo 



0(a' 



where 



Ai) 



Vk 

Ai) 



Ai) 

(i) 

W 

Ai) 



A/'(0,S(*)), 



1,2. 



(50) 



Correlation matrices B^^^ and B^"^^ in both sequences are strictly positively definite 
matrices. Denote by 71 and 72 the Lyapunov exponents for the sequences Mj^^^ and 



M 



(2) 



k = 1,2,.... Let det Mo 7^ 0, < a ^ ao, and ^ B^'^\ B'^^^ ^ A2/ with 
Ai^2 > 0. Then there exist constants Ci,C2 depending on Mo, \i, \2, and (Tq such that 
Ci72 ^ 7i ^ C272 uniformly with respect to Mq, B^^''^\ and a. 

The lemma has obvious physical meaning: the stronger the disorder in the waveg- 
uide, the greater is its Lyapunov exponent. However, mathematical proof of the lemma 
is quite tedious and will be published elsewhere. 
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5 The white noise model 



In the previous section we reduced calculation of the Lyapunov exponent to the product 
of random matrices 

M{u) = Mo[l + aVk + 0(^2)] , (51) 

where 

cos uj£ 



sinaj£ 



Mo= - , (52) 

—00 sm 001 cos 001 

and Vfc, k = 1,2,..., are independent identically distributed Gaussian matrices. The 
exact calculation of the Lyapunov exponent for the product of random matrices is, 
unfortunately, quite a challenging problem. However, it can be calculated for the 
Schrodinger equation with the white noise potential using the Ito calculus. To this 
end, we consider the model equation 

-ip'^x) + (Tw{x)'ip = oo'^ipix), t e R, (53) 

where w{x) is the standard white noise. Since w{x) is a generalized function, we have 
to define solution of the last equation. We approximate w{x) by a piecewise- const ant 
process 



X e [nA,{n + 1)A], n = 0,±l,±2. 



(54) 



where are independent normal random variables with = 0, Var^ = 1. The 
transfer matrix over the interval nA has the form 



n-l 



M{nA) = Yl 



k=0 



n-l 

n 

k=0 



COS 



^A 



oo sm 



A 



00^ 



cj2 _ £ife sin 



00^ 



sik, ^ 



cos 



00 



00^ 



1 ujA 

-00 A H 1 

oo 



O (A3/2) . 



(55) 



(56) 



If we introduce the polar coordinates 

^p{kA) = r{kA)sme{kA), 

ip'ikA) = ur{kA)cose{kA), 

then the values of ip {{k + 1)A) and oo^^tlj'{{k + 1) A) are related through the transfer 
matrix (55) as follows 

1 uA 



r((/e + l)A) 



sin^((A; + 1)A) 
cos^((A; + 1)A) 



-ooA H 1 



00 



r{kA) 



sme{kA) 
cos6'(A;A) 



r{kA) 



cos 9{k A) 



sin e{kA) + ooA cos 9{kA) 

CT^kVA 



-ooA 



oo 



sine{kA) 



(57) 
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Denoting 

we derive from (57) 

z{{k + 1)A) 



z{kA) = cot^(A;A) 



(5^ 



z{kA) - cuA 



a^kVA 



l + ujAz{kA) 

cr^kVA 



z{kA)-cuA 



[l-ujAz{kA) + 0{A^)\ 



{kA) -uA + - coAz\kA) + 0(A3/2). 

UJ 



(59) 



Taking in (59) the limit A — )■ we obtain a stochastic differential equation for the 
limiting phase 9{x) corresponding to the white noise potential 



dz{x) 



■ w{x) — UJ (1 + z'^{x)). 



dx UJ 

Similarly, for the amplitude r{x) we obtain from (57) 



(60) 



r^i{k + l)A) =r\kA) 



+ i cos 9{k A) + (^-uA 



r\kA) 



{sm9{kA) + ujAcose{kA)y 



sm9{kA) 
uj J j 



V 



1 + sin e{kA) COS e{kA) + sin^ d{kA) + 0( A^/^) 

UJ UJ^ 



(61) 



Taking the logarithm of both sides of (61) and sending A — )■ we obtain a stochastic 
differential equation for the amplitude r{x) 



dhirix) 



a 



sin 26{x) w{x) sin 6{x) cos 26(x) 



dx 2uj \ y \ y 2[j2 
Equation (60) defines generator £ of the diffusion process z{x) 



(62) 



UJ{1 + Z^) 



d 



2uj'^ dz"^ ^ ' ' dz 



(63) 



The distribution density P{t, zi, Z2) of the diffusion process satisfies the equation 

dP 



dt 



CP with „ = 5{zi - Z2). 



(64) 



The limiting distribution density p is then given by 

C*p = Q 



(65) 
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or 



2ijj'^ dz"^ dz 
Solution of this equation has the form 

p[z) = Ce"*(^) r e*« dt, $( 

J — oo 

where C is a normahzing constant 



M^ + y)^ ^ = ^^ (67) 



1 = r e-*(^) r e*W cit = dx. (6^ 

C y_oo J-oo V A Jo 



Now we can calculate the Lyapunov exponent 



liin Eln(r(L) ^ , jj^2^ 20d^ 



^2 /-oo ^_^2 



C / , e-'^(")dx / e*(^)dy. (69) 



2^2 (1 + X 

Rigorous mathematical analysis of dependence of 7 on a; and a is quite cumbersome 
and will be published elsewhere. Below we give the summary of the result. 



5.1 Frequency near the band edge 

Let (Jo 7^ be a band edge that separates a gap and a non-degenerate band (i.e. point 
A in Fig. 3), and let frequency u be in the band and close to the band edge uq such 
that leu — cjol = i^A ^ 1- Then for small a and 

2 3 r, 

7 = Cias provided < a , (70) 

where Ci is a constant. We illustrate this dependence numerically when a small disorder 
is introduced in the lengths li and I2 of the waveguide in Section 3. For propagation 
frequency u = 5.6288 corresponding to the left edge of the band (point A in Fig. 
3) we calculated the transmission coefficient after 1000 periods and then using (11) 
we obtained the Lyapunov exponent 7 averaging the last 500 periods. Results of the 
numerical simulation are shown in Fig. 4 by circles. The approximating line has 
equation 

lg7 = ^lga + 0.15 (71) 

which is in excellent agreement with the asymptotics (70). The same dependence was 
confirmed numerically in [7]. 



5.2 Frequency in the bulk of the band 

If frequency u is located in the band far from the band edges then the character of 
asymptotics of 7 is completely different. In this case 



7 = C2(T^ provided cr <^ 1. 



(72) 
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Figure 3: The discriminant curve and bandgap structure of model 1 with rii = 1, n2 = 
2.5, h = 1, I2 = 0.1. Point A corresponds to the left band edge u = 5.6288 while point 
B with u = 9 lies in the middle of the band. 

Here C2 is a constant independent of a. We illustrate this dependence numerically 
for frequency u = 9 (point B in Fig. 3). The result is shown in Fig. 4 by squares. 
Dependence of Ig 7 on Ig a is then approximated by the line 

lg7 = 21ga+ 1.34 (73) 

which is well consistent with the asymptotics (72). 

6 Conclusions 

We have studied the influence of disorder on transmission through periodic waveguides 
described by the Schrodinger and optical equations. While the bands of the Schrodinger 
operator have regular structure, it is shown that the bands of the optical operator do 
not possess this property, and the lengths of its gaps exhibits a chaotic behavior as the 
frequency Jncreases. Assuming that the transfer matrix of perturbed waveguide has 
the form Mk = M {I + aVk + 0{a'^)) we have shown that this representation is valid 
for two particular models of disordered waveguides. Using the result that Lyapunov 
exponents of systems with transfer matrices in the canonical form estimate each other, 
we have estimated the Lyapunov exponent of the optical waveguide through that of 
the Schrodinger equation with the white noise potential. Exact dependence of the 
Lyapunov exponent 7 was established as a function of frequency u and intensity of 
the disorder a. When the frequency u lies in the bulk of the band then 7 ~ o"^, while 
near the band edge 7 has the order 7 ~ cr^/^. Thus, small disorder drastically reduces 
transmission of the waveguide if the frequencies are located near non-degenerate band 
edges. 
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Figure 4: Numerical simulation of dependence of the Lyapunov exponent 7 on the 
strength of disorder a. Circles and squares correspond to perturbation of points A and 
B in Fig. 3, respectively, and approximated by the lines with the slopes 2/3 and 2. 
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